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ON  THE  STABILITY  OF  INTERACTING  QUEUES  IN  A  MULTIPLE  ACCESS  SYSTEM 

by 

R.  Rao 

University  of  California,  San  Diego 
and 

A.  Ephremides 
University  of  Maryland 
Abstract 

We  consider  the  standard  discrete-time  slotted  ALOHA  system  with  a  finite  number  of 
buffered  terminals.  The  stability  (ergodicity)  region  for  this  system  is  known  for  the  case  of  two 
terminals  and  for  the  case  of  infinite,  but  symmetric,  terminals.  In  this  paper  we  introduce 
a  new  approach  of  studying  the  stability  of  this  system  by  meanB  of  a  simple  concept  of 
dominance.  As  a  result  we  show  that  the  stability  region  for  the  case  of  two  terminals  can  be 
obtained  in  a  very  simple  way.  Furthermore,  we  obtain  lower  (inner)  bounds  for  the  stability 
region  of  the  system  with  an  arbitrary  finite  number  of  terminals.  These  bounds  are  superior 
to  the  ones  already  known.  Finally  we  point  out  a  similarity  betwen  these  stability  results 
and  the  achievable  region  of  the  no-feedback  collision  channel  that  may  suggest  a  connection 
between  the  two  problems. 

1.  Introduction 

Interacting  queueing  systems  occur  naturally  in  multiple  access  channel  models  and  shared 
computer  processor  systems.  As  is  the  case  with  most  non-standard  queueing  systems,  these 
interacting  queues  are  difficult  to  analyze.  Their  study  has  received  much  attention  lately 
owing  to  their  importance  in  applications,  as  well  as  to  their  theoretical  interest.  In  [1]  Fayolle 
and  Iasnogordski  displayed  the  inherent  difficulty  of  the  analysis  of  such  systems.  In  [2]  the 
importance  of  interfering  queues  in  multiple  access  systems  was  recognized.  In  [3]  Saadawi  and 
Ephremides  introduced  an  approximate  model  for  the  analysis  of  the  slotted  ALOHA  system. 
In  [4,5]  Sidi  and  Segall  introduced  different  approximations  that  led  to  an  exact  analysis  of  a 
simple  2-user  system.  Szpankowski  in  [6,7]  has  considered  the  ergodicity  region  of  the  slotted 
ALOHA  system  and  obtained  lower  bounds.  Tsybakov  and  Mikhailov  [8]  obtained  sufficient 
or  necessary  conditions  for  the  stability  of  this  system  for  the  case  of  M  users  and,  based 
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on  &  result  by  Malyshev  [9],  were  able  to  obtain  both  necessary  and  sufficient  conditions  for 
ergodicity.  Finally,  in  a  number  of  papers  that  have  appeared  in  the  literature  since  the 
early  seventies  it  has  been  established  that  in  the  classical,  bufferless,  symmetric,  infinite-user 
ALOHA  system  the  stability  threshold  is  t~l  [10,11,12]. 

In  this  paper  we  consider  the  case  of  discrete-time,  slotted  ALOHA  system  with  M  users 
each  of  which  has  a  buffer  of  infinite  capacity  to  store  incoming  packets.  The  assumptions  of 
the  model  are  the  usual  ones.  Time  is  slotted  and  the  transmission  time  of  a  packet  is  equal 
to  one  slot.  Each  user  receives  (generates)  packets  according  to  a  Bernoulli  process.  The  rate 
of  arrivals  is  A<  for  the  ith  user.  Arrivals  to  different  users  are  independent.  In  each  slot  user  i 
attempts  to  transmit  the  head-of-the-line  packet  with  probability  pi,  provided  the  buffer  is  not 
empty.  Based  on  instantaneous  ternary  feedback  (collison,  idle,  success)  each  user  determines 
the  outcome  of  the  attempted  transmission.  Simultaneous  transmission  attempts  by  two  or 
more  users  result  in  collision. 

We  are  interested  in  determining  the  region  of  values  of  the  arrival  rates  A<,«  =  1  ,...,M 
for  which  this  system  of  M  queues  is  stable.*  We  introduce  a  new  approach  in  the  study  of 
the  stability  question.  This  approach  consists  of  considering  hypothetical,  auxiliary  systems 
of  queues  that  closely  parallel  the  operation  of  the  system  of  interest  but  dominate  it  in  a  well- 
defined  sense.  As  a  result,  we  are  able  to  determine  the  stability  region  for  the  case  of  M=2, 
that  was  obtained  in  [8],  in  a  very  simple  way  that  illuminates  the  relationship  and  interaction 
between  the  two  queues.  This  is  done  in  Section  2.  Furthermore,  by  using  the  same  approach 
in  a  slightly  more  elaborate  form  we  are  able  to  obtain  improved  sufficient  conditions  for  the 
stability  of  the  system  of  M  queues,  for  any  finite  M.  This  is  done  in  stages.  In  Section  3 
a  first  set  of  bounds  is  obtained  from  the  direct  extension  of  the  dominance  concept  as  used 
in  Section  2  and  in  section  4  the  more  elaborate  set  of  dominating  systems  is  used  to  lead 

to  the  derivation  of  the  final  set  of  bounds.  Finally  in  section  5  we  relate  these  results  to 

0 

the  achievable  region  of  the  collision  channel  without  feedback  [13,14].  Some  of  the  technical 
proofs  have  been  moved  to  appendices  for  the  sake  of  clarity  and  continuity  in  the  main  text. 

2.  The  Dominant  System  for  the  Case  of  2  Queues 

Let  us  denote  by  S  the  system  of  queues  that  is  the  object  of  our  study  as  described  in 

*  the  precise  definition  of  stability  is  provided  in  the  next  section 
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the  introduction.  Let  Qi  denote  the  queue  size  at  the  ith  terminal,  t  =  1, 2.  Clearly,  the  queue 
sizes  evolve  as  a  2-dimensional  random  walk  in  the  first  quadrant  as  shown  in  Fig.  1.  The 
transition  probabilities  of  the  underlying  Markov  chains  can  be  easily  obtained,  based  on  the 
rules  of  operation  of  the  ALOHA  system.  In  fact  this  chain  is  irreducible  and  aperiodic.  What 
makes  the  analysis  difficult  is  of  course  the  infinite  size  of  the  state  space  and  the  fact  that 
the  probabilities  of  transition  on  the  boundaries  are  different  from  those  of  their  counterparts 
in  the  interior  of  the  quadrant.  This  is  shown,  for  example,  in  the  figure  for  the  transition 
(Q11Q2)  — ►  (Q 1  — l,<?a)-  If  these  transition  probabilities  were  identical,  regardless  of  whether 
Q2  =  0,  the  system  of  two  queues  would  decompose  nicely*  and  each  queue  could  be  separately 
analysed.  This  observation  motivates  the  introduction  of  an  auxiliary  hypothetical  system  that 
has  precisely  this  property.  Consider  a  new  system,  S2,  consisting  of  copies  of  the  two  queues 
of  5,  with  the  following  properties: 

1)  arrivals  at  queue  i  in  the  new  system  occur  at  exactly  the  same  instants  as  in  the 
original  system,  <  =  1, 2; 

2)  the  "coin-tosses”  that  determine  transmission  attempts  at  queue  »,»  =  1,2,  have 
exactly  the  same  outcomes  in  both  systems; 

3)  whenever  Qi  =  0,*  =  1,2,  terminal  »  continues  to  transmit  “dummy”  packets  with 
the  same  probability  p*,  thus  continuing  its  interference  with  the  other  terminal, 
whether  it  is  empty  or  not. 

It  is  clear  that  the  queue  sizes,  at  either  terminal,  in  the  new  system  will  never  be  less  than 
their  counterparts  in  the  original  system,  provided  they  start  from  identical  initial  conditions 
in  both  systems.  In  the  new  system  5 2  the  “service”  rate  seen  by,  say,  terminal  1  is  always 
equal  to  pj(l  —  p2)  while  in  the  original  system  it  oscillates  between  the  values  of  p%  and 
Pi(l  —  p2)  depending  on  whether  terminal  2  is  empty  or  not. 

Definition.  The  stability  region  of  either  system  is  the  set  of  values  of  A,-,  t  =  1,  %Sot  which 
the  underlying  irreducible,  aperiodic  Markov  chain  is  ergodic. 

Clearly,  the  strong  stochastic  dominance  of  S2  on  S  implies  that  the  stability  region  of 
S2  inner-bounds  that  of  S.  The  stability  region  of  52  is  easy  to  determine.  Each  terminal 

*  Note,  however,  that  the  queues  would  still  be  statistically  dependent  since  the  transition 
[QuQz)  — ►  (Qi  —  1,Q2  —  1)  has  zero  probability. 
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operates  as  a  discrete  M  |  M  |  1  system  and  is  therefore  stable  if  and  only  if 

xi  <  Pi  11(1 " ») 
jY* 

Thus  the  stability  region  of  S3  is  the  rectangle  shown  in  Fig.  2.  This  is  of  course  a  poor  bound 
for  the  stability  region  of  S  because  53  dominates  the  original  system  too  strongly.  We  would 
like,  to  draw  attention  to  one  observation  before  proceeding  any  further.  Note  that  each  queue 
behaves  identically  in  both  systems  so  long  as  neither  queue  empties.  Thus  if  both  systems 
get  started  with  identical,  non-zero  queue  sizes,  they  will  evolve  following  exactly  identical 
trajectories  at  least  until  one  of  the  queues  empties.  This  observation  plays  a  crucial  role  in 
some  of  the  proof  arguments  that  will  follow. 

Now,  let  us  consider  another  hypothetical,  auxiliary  system  S1,  in  the  hope  of  obtaining 
a  dominant  system  that  does  not  dominate  the  original  system  as  strongly  and  that  can  lead, 
therefore,  to  a  closer  tracking  of  the  original  system’s  behavior  and  to  a  better  bound  of  its 
ergodicity  region.  The  system  S 1  is  identical  to  both  S  and  5s  with  respect  to  arrivals  and 
attempted  transmissions  (that  is,  it  has  properties  (l)  and  (2))  but  differs  with  respect  to 
the  behavior  when  the  queues  empty.  In  this  system  queue  2  behaves  as  in  S  while  queue  1 
behaves  as  in  S3;  that  is,  only  queue  1  continues  the  irresponsible  transmission  of  “dummy” 
packets  when  it  empties.  As  a  result,  queue  1  sees  a  “service”  rate  that  oscillates  between  pi 
and  pi(l  —  Pi)  depending  on  whether  queue  2  is  empty  or  not,  while  queue  2  always  sees  a 
worst-case  service  rate  of  pj(l  —  px)  regardless  of  the  status  of  queue  1.  It  is  clear  that  S1 
dominates  S  since  either  queue  will  have  a  successful  departure  in  S  whenever  it  has  one  in 
S1,  but  not  necessarily  vice-versa.  It  is  also  clear  that  in  Sl  the  stability  of  queue  2  is  easily 
determined  by  the  fact  that  it  operates  as  a  discrete-time  M  \  M  \  1  system.  Thus  it  will  be 
stable  iff 


<  Ps(l  —  Pi)  (1) 

Let  us  therefore  assume  that  (l)  is  satisfied  and  let  us  determine  the  criterion  for  the 
stability  of  queue  1.  Terminal  1  sees  a  “service”  rate  that  has  the  value  pi  when  terminal  2  is 
empty  (and  that  happens  with  probability  1  -  according  to  the  well-known  M  |  M  |  1 
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formula)  and  the  value  pi(l -p2)  when  terminal  2  is  non-empty  (which  occurs  with  probability 
pTC^-Pil).  ^PP611^**  ^  “  shown  that  a  necessary  and  sufficient  condition  for  the  stability 
of  queue  1  is  that  the  arrival  rate  Ai  be  less  than  its  average  service  rate,  thus  extending  the 
well-known  result  for  a  G  |  G  |  1  queue  [15].  Applying  this  result  here  yields 

Al  <  *(1  ’  + ',l(1  - 
which  can  be  written  as 

A,  <  ,,(1  -  -iM  (2) 

1  -Pi 

Thus  the  stability  region  of  S1  is  as  shown  in  Figure  3. 

What  is  important  is  that  the  boundary  of  the  stability  region  given  by  Eq.  (2)  is  not 
only  a  bound  for  the  stability  region  of  S  but  coincides  with  it.  The  proof  relies  on  the 
observation  made  earlier  that,  so  long  as  the  queues  do  not  empty,  systems  S  and  S1  are 
indistinguishable.  Here  is  the  argument:  Given  that  A2  <  p3(l  —  pi),  if  for  some  A*  queue  1  is 
stable  in  S1,  it  is  also  stable  in  S  by  virtue  of  the  dominance.  Conversely,  if  for  some  Aj  queue 
1  is  unstable  in  S1,  then  the  queue  size  Qi(n)  will  grow  to  oo  without  emptying  with  finite 
(non-zero)  probability.*  Thus,  not  all  sample  paths  of  Qi(n)  correspond  to  transient  behavior 
with  infinitely  often  visits  to  0.  However,  so  long  as  queue  1  does  not  empty,  51  and  S  behave 
identically  if  started  from  the  same  initial  conditions.  Thus  the  same  sample  paths  that  go 
to  oo  without  visiting  0  belong  to  the  evolution  of  queue  1  in  system  S.  Therefore  queue  1  is 
unstable  in  system  5  as  well. 

Notice  that  by  reversing  the  roles  of  the  two  queues  in  system  51,  that  is  by  assuming 
that  queue  2  is  the  one  that  transmits  dummy  packets  when  it  empties,  we  obtain  &  stability 
region  given  by 

Ai  <  pi(l  -  pa) 

*  This  follows  from  the  fact  that  queue  2  is  stable  in  its  own  right  and  decoupled  from  queue 
1;  nevertheless  an  independent  proof  can  be  supplied  as  in  [16]. 
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and 


A,  <  P,(l  -  —-)  (3) 

1~P2 

The  branch  of  the  boundary  that  corresponds  to  Eq.  (3)  can  be  shown  to  be  part  of  the 
stability  region  of  5  in  a  similar  fashion.  Thus  the  union  of  the  branches  corresponding  to 
Eqs.  (2)  and  (3)  defines  the  stability  region  for  S  as  shown  in  Fig.  4.  This  is  identical  to  the 

region  obtained  by  Tsybakov  and  Mikhailov  in  [8]  through  use  of  Malyshev’s  lemma.  Note 
that  by  taking  the  envelope  of  these  regions  as  pi  and  pi  vary  in  [0, 1]  we  obtain  the  curve  C, 
shown  in  Fig.  4,  which  is  analytically  described  by 


y/\x  4-  \Z^a  =  1 

or  parametrically,  by 


(4) 


Ai  =  Pi(l-Pa) 


^2  =  Pa(l  —  Pi) 


(5) 


Pi  +  P2  =  1 

The  idea  of  this  simple  derivation  as  well  as  the  choice  of  dominant  systems  was  first 
briefly  reported  in  [16]. 

3.  The  Dominant  System  for  the  case  of  M  queues 

Now  consider  the  case  of  M  >  2.  The  approach  described  in  section  2  is  insufficient  to 
yield  the  ergodicity  region  of  the  system.  It  does  help  however  in  producing  bounds  that  are 
superior  to  the  ones  known  to  date.  In  this  section  we  will  develop  this  dominance  approach 
and  show  how  these  bounds  can  be  obtained. 

Let 

Qi(n)  =  queue  size  at  terminal  i  at  time  n 
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2?i(n)  =  number  of  departures  from  terminal  i  at  time  n,»  =  1  ,...,A/,n  =  1,2,.... 


and  let 


p$ = p?m«) = *] 


M»(n)  =  Pr[Di(n)  =  l]  (7) 

Thus  /*i(n)  is  the  “service”  rate  at  terminal  i  at  time  n.  It  clearly  depends  on  the  status  of 
the  queues  at  the  different  terminals  and  assumes  values  between  p»,  when  all  other  queues 
are  empty  and  p,  11^(1  ~  Pj)>  when  none  of  the  other  queues  are  empty.  Let  £*(n)  denote 
the  average  value  of  m(n).  We  further  define  two  types  of  events  of  interest: 

£*(n)  =  {<2*(n)  =  0},  i.e.  the  event  that  at  time  n  the  kth  queue  is  empty  (8) 


and 


Ekl *,. „*.(*)=  n  n  &W  (®) 

i.e.  the  event  that  at  time  n  the  only  empty  queues  are  the  ones  corresponding  to  termi¬ 
nals  fce.  Also  note  that  the  family  of  everts  E  forms  a  partition  of  the  sample  space. 

Furthermore  we  have 


&(»)  =  U  ^**»(w)}U-LK  U  Ekk' (»)) U U 

*i#* 


*i . *. 

•11  dlitiact 


(10) 


In  Appendix  1  it  is  established  that  the  limits  of  m»(«)  as  »  -»  oo  exist;  we  denote  them 
by  fix  and  their  average  values  by  fo.  Let  us  consider  also  the  system  S1,  defined  in  a  manner 
similar  to  that  of  the  preeeeding  section.  Namely,  S'1  is  identical  to  S  in  all  respects  (arrival 
instants,  transmission  attempts,  etc.)  except  that  terminal  1  continues  to  transmit  “dummy” 
packets,  even  when  it  is  empty,  with  the  same  probability  pi  as  when  it  is  not  empty.  All  other 
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terminals  behave  as  in  5.  We  define  the  quantities  Q1 1 D1 ,  P1 ,  and  fi1  as  above  but  with  the 
use  of  the  superscript  "1”  to  denote  reference  to  system  S1.  Similarly  we  also  consider  for 
reference  purposes  the  system  SM ,  in  which  all  terminals  continue  to  transmit  when  empty. 
Thu  system  is  the  generalization  of  S*  as  defined  in  Section  2. 

Our  goal  is  to  derive  sufficient  conditions  for  the  stabilityof  Sl  and  to  conclude,  by  virtue 
of  the  dominance  relationship,  that  these  are  also  sufficient  conditions  for  the  stability  of  S. 
We  shall  focus  on  queue  1.  To  "decouple”  its  stability  status  from  that  of  the  other  queues 
we  need  to  assume  that  the  other  queues  are  stable.  Unfortunately,  unlike  the  case  of  M  —  2, 
we  do  not  know  what  is  the  stability  criterion  for  the  set  of  queues  2, ...,  M  and  thus  we  are 
forced  to  assume  a  "strong”  stability  status  for  them,  namely  that 

<  Vi1  =  Pi  JJC1  ~  Pj)>  *  =  2, M  (11) 

The  validity  of  Eq.  (10)  implies  stability  of  this  group  of  queues  under  system  SM  and  thus 
under  S1  and  S  as  well. 

Let  B\(n)  denote  the  event  that  in  system  S1,  terminal  1  at  time  n  successfully  transmits 
a  packet,  whether  real  of  “dummy”;  its  probability  is  /i](n).  We  shall  presently  show  that  a 
sufficient  condition  for  the  stability  of  queue  2  in  S1  is 

Ai  <  ^(Aj,...,  Am; Pi,—, Pm)  (12) 

where  ^  is  a  specific  function  of  the  other  arrival  rates  and  all  retransmission  probabilities, 
that  lower-bounds  the  asymptotic  service  rate  Mi  of  terminal  1.  In  Appendix  2  we  show  that 
this  bound  is  already  an  improvement  over  these  derived  in  [6-8].  We  present  the  results  in  the 
form  of  two  lemmas  that  lead  to  a  theorem.  The  first  lemma  lower- bounds  ^J(n),  the  second 
lemma  establishes  the  bound  as  n  — ►  oo,  and  the  theorem  establishes  the  sufficient  condition 
by  making  use  of  the  results  of  Appendix  1. 

Lemma  S.l 

In  the  system  S1,  the  probability  of  a  successful  transmission  by  queue  1  at  the  instant  n 
satisfies  the  following  lower  bound  inequality  for  all  »  >  1: 
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(13) 


M  M  M 

tA  n  >  pi  n  “ ») + Pi  a  pr\ $  (n)i  •  Pi  n  u  ~ ») 


j=2 


3=2 


y«=s 

i#. 


Proof 

By  definition: 


„!H  =  Pr[B‘(n)|  (14) 

We  apply  the  theorem  of  total  probability  to  Eq.  (14)  conditioning  the  event  J3j(»)  on  the 
2m~1  mutually  exclusive  events  of  Eq.  (9).  This  results  in  the  following  equation: 


#•!(")  =  Pr[B\ (n)]  =  Pr[Bj(n)/{^(n)]  •  Pr [«}(»)]+ 

+  £;  Pr[B\  ( n)/Bl,  (»)]  •  Pr|El,  (»)]  + ... 

*i=3 


...+  £  Pr[B!  (»)/<...,,(«» -PflE*. . *.(«)]  + 

»••••  *««{2  »•••» 
mil  kj  diatimct 


+Pr[B{  (n)  |  Ba(...^#(n)]P»’[Ba,...>/{n)] 


(15) 


The  following  further  decomposition  is  now  possible. 

Pr[Bj(»)/Bj . (»)] 

=  PrIBjW/WiW  >  0 },El . „»]  •  Pr[{<3l(n)  >  0}/P] . ,.(»)]+ 

+Pr[B,' (»)/{<?!(»)  =  0).  El, . ,.(»)]  ■  Pr|(«}(n)  =  0>/Ej, . *»)  (16) 


B 


Now,  if  queue  1  is  non-empty  and  of  the  remaining  ones  only  queues  h\  through  kt  are  empty 
then  the  probability  of  a  successful  transmission  by  terminal  1  is  p\  n»^{*lt...pk.}(l  ~  P«)-  Note 
that  in  this  case  the  transmitted  packet  is  a  real  information  packet.  On  the  other  hand, 
if  queue  1  is  empty  and  the  status  of  the  other  queues  is  as  before  then  the  probability  of 

a  successful  transmission  is  still  given  by  pi  II*>f  . *,}(*  —  P»)-  The  difference  is  that  in 

this  case  the  transmitted  packet  is  a  dummy  packet.  However,  we  are  not  concerned  yet 
with  this  distinction  since  at  this  stage  we  seek  a  lower  bound  for  the  probability  of  successful 
transmission  of  information  or  dummy  packets  n\(n).  Besides,  near  saturation,  all  packets 
tend  to  be  information  packets  since,  then,  the  queues  are  rarely  empty. 

Based  on  the  above  comments  we  deduce  that 

Pt[B\  («)/£;„...,.(»)]=  Pi  n  (» -  P.)  (W) 

Now,  note  that  the  following  inequalities  can  be  easily  verified  for  any  set  of  indices 
fcj,  fcj,  ,..kt  i 

1  =  n  G>* + (i - »)) = (p*» +i-p*»)  n  b* +(!-?>*))= 


=[«. +(!-«,)]  n  u-!>*)(i+rr^)>»,  n  (»-»)+ 

**{*»,*»  }  *<{*a.*» 


+  n  (»-»)  ..  <is) 


i = n>* + c1  ~  p*)i  ^  m1-*) + n  »» 

k=2  k=2  *i=2 


(19) 


Using  eqns.  (17),  (18)  and  (19),  we  rewrite  Eq.  (15)  as  follows: 
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M 

/*Kn)  >  pi  n  ~  w)  •  pr[£}(n)i+ 


*= 2 


+  53  k*  no- p*)) + {pip*i  no- **))]  •  •pf[p*1  w) + — 


*i=2  *= 2 


k=2 

»**l 


...  + 


M 


M 


i{pi  n  o  -  **)> + (pip*i  no-  p0)>]  •  p,,[p*i  ,*„...*.  w] + - 


<k1.ka....k,)t<3.#,...kO  *=2 

all  kj  diitiict 


ks3 


j if 

- + kpi  no- **)> + *  53  pip*i  no-  p*)}]p*-[£2\s,...Af  w] 


*= 2 


*1=2 


k=3 

M“l 


(20) 


Rearranging  Eq.  (20)  and  using  eqs.  (9)  and  (10)  we  obtain  the  following  inequality 


m!(«)  >  Pi  n0  “Pi)  +  53*{*  HO  -Pi)}Pf[f‘(n)]  (21) 

i=2  i=2  g> 

This  concludes  the  derivation. 

In  the  next  result,  we  derive  the  asymptotic  value  of  the  bound  of  Lemma  3.1.  The  bound 
b  calculated  on  the  assumption  that  queues  2  through  M  are  stable. 

Lemma  3.2 
If 


A*  <  Mi*  =  Pi  HO  “  Pi)» 
jV* 


»=1 . M 


then 


J™,  (n)  -  *  n  0  -  Pi)  +  53  **0  -  -nt? )  n  0  ~  Pi) 

“•  y=3 


i=2  »=2 
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Proof: 

It  is  sufficient  to  show  that 

lim  Pr[£*  (n)]  >1  — ^  2  <  *  <  M 

n-»oo  fl?*  ““ 

Note  that  by  definition 

^(n)]  =  Pr[Qj(n)  =  0] 

Since  SM  dominates  S1  we  know  that 

P|«:(")  =  0]>Pr[g“(n)  =  0] 

Therefore,  taking  limits  and  using  the  decoupling  properties  of  SM  we  obtain 

Jim  Pr\C; (n)]  >  Jim  Pr[Q,M(n)  =  0]  =  (1  -  A)  (22) 

Taking  limits  on  eq.  (21)  and  employing  the  inequality  (22)  we  obtain  the  desired  result, 
namely 


J!®,  A»i(n)  ^  *  II  (X  “  Pi)  +  Pi  X^1  ~  T^rK1  ~  Pi)  =  A  (23) 

3=2  *=2  P* 

The  existence  of  these  limits  is  again  justified  in  Appendix  1. 

Now  we  proceed  to  establish  the  main  theorem. 

Theorem  3.1 

The  system  S  is  stable  if 

M  M  hi 

<  p.  ri(i  -  ?,) + pi  ep,(i  -  4)  n  (‘  -  p*)  p4) 

i=J  t- J  “i  S’ 

and 
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2  <  j  <  M 

where 


M 


=pj  n(*-»> 

h=i 


Proof 

Prom  Lemma  3.2  we  know  that  the  R.H.S.  of  Eq.  (24)  is  a  lower  bound  for  p*,  hence 
Ax  <  fi{  and  Ax  <  jx\.  According  to  the  results  in  Appendix  1,  this  is  sufficient  to  ensure  that 

Pr[Q\  =  0]  >  0 

The  inequalities  (25)  are  also  sufficient  to  imply  that 

Pr[<#  =  0]  >  0,  2  <%<M 

Thus  under  the  hypothesis  of  the  theorem  the  system  S1  is  stable  and  so  is  5  since  it  is 
dominated  by  S1. 

Comments 

1.  Recall  that  in  System  S1  we  required  that  terminal  1  transmit  dummy  packets  when 
empty;  the  choice  of  terminal  was  clearly  arbitrary.  If  we  interchange  the  roles  of 
terminal  1  and  terminal  i  and  repeat  the  proof  we  will  derive  a  similar  result  for  a 
different  ordering  of  the  terminals.  In  fact  the  union  of  the  regions  that  are  obtained 
in  this  manner  for  all  possible  orderings  produces  a  better  bound. 

2.  Szpankowski’s  results  on  the  stability  of  system  S  [6,7]  are  of  similar  nature.  In 
appendix  2  we  show  that  the  stability  region  established  in  [6,7]  is  a  strict  subset 
of  the  region  obtained  here.  We  should  like  to  record  that  Szpankowski’s  results 
apply  to  systems  for  which  the  arrival  rates  are  not  necessarily  Bernoulli,  and  in 
that  sense  his  results  are  more  general.  Nevertheless,  we  believe  that  our  results 
can  also  be  easily  generalized  to  accomodate  non-Bernoulli  arrival  processes  because 
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the  dominating  systems  S1  and  SMt  crucial  to  our  analysis,  retain  their  essential 
properties  for  arbitrary  arrival  processes. 

3.  The  bound  obtained  in  Theorem  3.1  is  not  satisfactory.  As  shown  in  Fig.  5,  it  fails 
to  provide  values  for  regions  of  values  of  the  A,’s  that  are  outside  the  strict  strong 
stability  ones  implied  by  Eq.  (25).  These  regions  correspond  to  the  areas  marked  by 
question  marks  in  Fig.  5.  Also  it  is  clearly  not  very  tight  since  on  the  A,-axes  the 
actual  stability  surface  intersects  at  the  points  pi,P2,p$  etc.  In  the  case  of  M  =  3, 
as  shown  in  Fig.  5,  the  bound  is  below  these  values  by  an  amount  equal  to  pipaps. 

4.  A  Series  of  Dominant  Systems  that  Yield  Improved  Bounds 

In  the  preceeding  section  we  obtained  an  inner  bound  to  the  stability  region  of  S  by 
essentially  using  the  same  approach  as  in  section  2  in  which  we  were  successful  in  determining 
exactly  the  stability  region  for  the  2-user  system.  In  this  section  we  refine  our  approach  by 
considering  a  series  of  dominant  systems  that  are  able  to  track  the  behavior  of  S  a  little  closer 
then  S 1  could.  These  systems  permit  us  to  by-pass  the  shortcomings  of  the  bound  of  Theorem 
3.1  and  to  relax  the  sufficient  conditions  for  the  stability  of  5. 

Retaining  the  notation  we  have  used  so  far  we  introduce  now  an  extended  set  of  systems. 
For  j  =  1, 2, M  we  define  S3  as  follows: 

1)  Arrivals  at  the  ith  queue  of  S3,  *  =  1, ...,  M,  are  identical  to  those  at  the  ith  queue  of 

5. 

2)  “Coin  tosses”  that  determine  transmission  attempts  at  the  ith  queue  of  S3,i  = 
1, ...,  M,  have  identical  outcomes  to  those  at  the  ith  queue  of  S. 

3)  For  i  >  j,  terminal  i  behaves  exactly  as  in  S,  that  b  it  does  not  attempt  to  transmit 
“dummy”  packets  when  empty. 

4)  For  *  <  j,  terminal  i  attempts  to  transmit  “dummy”  packets  when  empty  according 
to  the  following  rules:  with  the  aide  of  a  “genie”,  terminal  i  b  informed  whether 
any  terminal  k,  with  k  <  t,  will  attempt  a  transmission  in  the  slot;  if  yes,  terminal 
i  refrains  from  attempting  to  transmit;  if  no,  it  attempts  to  transmit  a  “dummy” 
packet  with  probability  p,-. 

As  a  result  we  have  a  series  of  systems  in  each  of  which  terminal  1  b  afforded  the  first 
choice  to  transmit  a  “dummy”  packet,  if  empty.  If  it  faib  to  transmit,  “dummy”  or  real 
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packet,  terminal  2,  if  empty,  is  afforded  the  same  choice.  If  it  doesn’t  transmit,  terminal  3 
is  given  the  chance,  and  so  on,  up  to  terminal  j  in  the  system  S’.  Terminals  with  identity 
number  greater  than  j  do  not  attempt  to  transmit  “dummy’’  packets  in  S’. 

As  a  result,  this  series  of  systems  has  the  following  properties. 

i)  S’  dominates  S’~ltj  =  2,...,  M,  and,  of  source,  S1  dominates  S.  Note  that  S 1  coin¬ 
cides  with  the  system  S 1  considered  in  section  3  and  so  does  SM  with  its  counterpart 
in  that  section.  The  reason  for  the  dominance  is  this:  Terminal  j  in  S’ ~l  is  always 
silent  when  empty,  while  it  may  transmit  a  “dummy’’  packet  in  S’.  Furthermore,  all 
other  terminals  in  S’~l  will  have  a  successful  transmission  whenever  they  do  have 
one  in  S’ ,  but  not  vice-versa.  Thus,  it  is  clear  that 

<#(n)  >  Qr'in)  >  Qi[n)i  »  =  1,2,...;/  =  2,...,M;»  =  1  ,...,M 

ii)  The  joint  queue  sues  of  the  M  users  in  each  of  these  systems  evolve  according  to  an 
M-dimensional,  aperiodic,  irreducible  Markov  chain. 

iii)  In  system  S’  the  cummulative  contention  from  terminals  1,2,...  j  for  the  use  of  the 
channel  remains  fixed  and  equal  to 


i-no-K) 

»=1 

irrespective  of  whether  any  of  these  queues  are  empty.  Thus  any  terminal  k,k  >  /, 
faces  competition  from  the  group  of  the  first  j  terminals  that  stays  always  fixed.  This 
is  ensured  by  the  rules  of  transmission  explained  above  in  (4). 

We  present  the  derivation  of  the  new  bounds  in  the  form  of  some  lemmas  that  lead  to  a 
theorem.  We  start  by  introducing  some  notation  similar  to  that  of  Eqs.  (8),  (9),  and  (10). 

For  any  set  of  distinct  indices  (/  +  1)  <  K  <  M ,  let 
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«>)  =  {«U»)=  o} 

<(»)={«;.(»)=“>  u  <«;(»)  > «} 

t=i+i 
t*  ki* 

k,m . *.(»)=  n  {«jw= 0}  n  m;m  >  °)  (26) 

•n*l**8*—*#} 

E}( »)  =  n  «?(»)>«> 

*•=;+! 

Thus  (n)  is  the  event  that  the  Jfcith  queue  of  SJ  is  empty  at  time  n.  E{  {n)  is  the  event 
that  only  the  Jkith  queue  of  SJ  in  the  group  from  j  +  1  to  M  is  empty  at  time  n  (all  others 
being  non-empty).  (n)  is  the  event  that  only  queues  fci, Jfc2, -■.£*  are  empty  at  time 

n  and  all  the  others  in  the  group  are  non-empty.  Since  the  group  of  terminals  from  1  to  j, 
collectively,  present  the  same  behavior  in  competing  for  the  channel,  regardless  of  the  status 
of  their  queues,  there  is  no  interest  in  whether  they  are  empty  or  not.  Thus  we  consider  only 
the  events  that  pertain  to  the  group  {/  +  1,  •  ••j  M). 

There  are  events  of  the  type  E3ki  (n)  and  (**”*)  events  of  the  type  E*ki  klt  kt{n). 

The  total  number  of  events  is  2M_J,  and  they  are  all  mutually  exclusive.  Their  union 
exhausts  the  sample  space.  The  event  £l  may  be  decomposed  into  mutually  exclusive  events 
as  follows. 


«U»)“«»)U<  U  <>>»•••  (27) 

*a  «{/+i....*f> 

*a**i 

We  now  prove  a  lemma  that  generalizes  lemma  3.1. 

Lemma  4.1 

In  the  system  SJ\  the  probability  of  a  successful  transmission  by  terminal  j  at  time  n 
satisfies  the  following  lower  bound  inequality  for  all  n  >  1, 
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Proof 


M 


i=l 


M  M 

E  *•[«/(»)]*  n(‘-«) 

•=i+ 1  k=y 


By  definition, 


/-;(»)  =  JV[B}(»)]  (28) 

where  Bj  (n)  ia  the  event  that  at  time  n  terminal  j  has  a  successful  transmission  of  a  real 
or  a  “dummy”  packet  in  system  j.  We  apply  the  theorem  of  total  probability  to  Eq.  (28) 
conditioning  the  event  B|(n)  on  the  2M~*  mutually  exclusive  events  of  Eq.  (26).  This  results 
in  the  following  equation. 


Fr[fl>  (»)]  =  Pr\Bi (»)/{£$(»)«  •  Fr[B>  („)]  +  £  P,[B>  (»)/£’,  (n))  -  Pr[<  („)]  + 


-P'K  Aij.(«)] + 


*< 


all  dUtlact 


+Pr[B>  (n)  |  E’+l . M(n)) .  Pr[£>+1 . M(n)]  (29) 

A  successful  transmission  by  terminal  j  is  impossible  if  any  of  the  stations  1  through  j-1 
transmit.  Let  us  define  by  Tf  (n)  the  event  that  terminal  i  attempts  a  transmission  in  slot  n 
in  system  j.  Then  we  decompose  the  general  term  of  Eq.  (29)  according  to  the  following: 


j-i 


PrlB’M/Ei^ (n) J  =  Pr[BJ(»)/  f|  WM}' flWiM  > 


t-1 


3-1 


■Pr !  n  {r/(»)}7W?(»)  >  0}  n  Ei . *.(”)]  •  >  <>/<,*„...*»)+ 


»=1 


17 


i-i 


+Pr\B‘(n)/  f|  (WF fKflJW  =  0}  H ».(«)]• 

3-1 


3-1 


■p'i n w(’*)}‘/{«?(") = o} n < . ».(»)]  •  ^K<?jw = . *.(")]  (30) 


t=i 


Since  the  competition  from  the  first  group  of  j-1  terminals  is  fixed  we  have 


nWwr/wjw = . *.w] = no  -  »>  oo 

»=i  k=i 

The  same  result  holds  for  the  case  of  Q^n)  >  0.  Now  let  us  consider  the  first  term  of  the  first 
summand  of  Eq.  (30).  The  conditioning  event  implies  that 

i)  terminals  1  through  j-1  remain  silent  at  time  n. 

ii)  terminal  j  has  a  non-empty  queue  and  hence  will  independently  attempt  transmission 
of  an  information  packet  with  probability  pj. 

iii)  of  the  remaining  terminals  only  Jfci  through  ke  are  empty  and  hence  they  will  remain 
silent  (according  to  the  rules  of  system  SJ). 

Based  on  this  we  deduce  that 


i-i 


?'{b’  (»>/  n  wm}'  nw;:M  > »}  n  <.» . »»i 


i=l 


=  Pi  II  (!-»)  (32) 

<«<y+» . u) 

Consider  now  the  first  term  of  the  second  summand  of  Eq.  (30).  The  conditioning  event 
implies  that 

i)  terminals  1  through  j-1  remain  silent  at  time  n. 
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ii)  terminal  j  has  an  empty  queue;  therefore,  by  the  rules  of  in  S*t  it  will  broadcast  a 
dummy  packet  with  probability  py. 

iii)  Of  the  remaining  terminals  only  Jbi  through  ke  are  empty  and  hence  they  will  remain 
silent.  Based  on  this  we  have 


nwwj'nwiw = o> o; . ».h]= p,  n  <> - w 

t=i  <=y+i . ** 

. . . 

Note  that  when  the  j*h  queue  is  empty  terminal  j  broadcasts  a  dummy  packet,  whereas 
when  the  queue  is  non-empty  it  broadcasts  an  information  packet.  However,  this  distinction  is 
again  of  no  concern  since  we  seek  a  lower  bound  for  the  probability  of  successful  transmission  of 
information  or  dummy  packets,  and  since  the  distinction  becomes  meaningless  near  saturation. 

Substituting  Eqs.  (31),  (32)  and  (33)  in  Eq.  (30)  we  obtain 

Pr[B’{n)/EiU  'k'(n)]  =  Pj  J]  (1  -  Pi)  (34) 

<*«>i . *.) 


Also  note  that 


Pr[BJ  («)/£'(»)]=,,,  J]  (1-p.)  (35) 

Note  that  the  special  nature  of  S3  enabled  us  to  obtain  eqns.  (34)  and  (35)  without 
having  to  estimate  Pr{Qj(n)  >  0 /£*,„..*.(»)]  or  Fr[Qj(n)  =  0/E^  .  ^(n)],  which  would  not 
be  possible. 

Now  we  invoke  the  inequality  (18)  and  a  slightly  modified  version  of  inequality  (19), 
namely, 


M 


M 


M 


M 


i=  n  ip* + (i - p*)i >  n  (i-p*)+  s  n  ^ 

*=j+ i  *=j+i 
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Substituting  eqns.  (34),  (35)  and  using  inequalities  (18)  and  (36)  in  eqn.  (29)  we  obtain 


JrL  mm  m 

My(n)>PiJ l(1-P*)-Pr[^(n)]+  £  {Pj  i1  ~P*)  +  PjPkx  fl  (1  —  (»)]+ 

S}  *i=j+i  *=».  *=>  ‘ 


W  W 

E  (foil  (»-»))  +  fop*.  II  (»-»)}) . ,.(»)]+ 


•  II  dUttect 


*=1 

Mi 


tel 

M*i.i 


W  •  M  M 

•••+({pJ  n(i-p*)>+{  e  *»;*>*.  n  u-p*)})-^^  M(„)j 

Mi  Al=>+1 


(37) 


or,  by  further  rearranging  and  using  Eq.  (27), 


M 


M 


M 


<‘j(n)aPiII(1-P')+  E  p)<Pi  5  <»-»)>*«(•)] 

*=i+i 


i=l 

Mi 


tel 

M<.i 


(38) 


This  concludes  the  derivation. 

By  letting  »  — ►  oo  in  Eq.  (38)  we  obtain  the  following  asymptotic  bound. 

M  M  M 

%  >  pj = pjw  >  pj  n  Pi  +  e  o»Pi  n  rjjmgj = «] 


:«  “<+> 


(39) 


k=l 

M»'.i 


It  is  clear  now  that  lower  bounds  for  n3-  can  be  found  if  lower  bounds  for  Pr[Q\  =  0]  can 
be  determined.  The  latter  are  derived  in  the  next  lemma  by  examination  of  the  departure 
statistics  at  queue  i  in  S3  for  j  <  *,  and  assuming  that  A i  <  p\t  for  %  >  /,  where  /2j  is  the 
average  asymptotic  service  rate  of  queue  i  in  system  S*. 

Lemma  4.2 

1£  Xi  <  for  all  i,  then  in  steady  state  Pr[Ql  >  0]  <  for  any  j  <  ».  * 

*  This  result  is  similar  to  the  one  for  the  G  |  G  |  1  queue,  see  Kleinrock  [15],  with  average 
arrival  rate  A  and  average  service  time  x.  If  the  queue  is  stable  for  A  <  |,  then  Pr[Q  >  0]  =  Ax. 
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Proof: 


Let  X?  (n)  denote  the  number  of  packets  in  the  ith  queue  of  S3  after  the  n*h  departure. 

Let  A\  (n+l)  denote  the  number  of  packets  that  arrive  at  queue  i  of  S3  during  the  (n+l)** 
service  interval.  Then: 


X’M  +  1)  =  X’(n)  +  A’(n  +  1)  -  ![*/(«)]  (40) 


where 


x  >  0 
else 


We  know  from  Appendix  1  that  if  A;  <  then  the  ith  queue  of  S3  is  stable.  Hence 
there  exists  a  steady  state  distribution  for  the  system  satistics.  Furthermore,  the  results  of 
Appendix  1  imply  that  the  average  queue  sizes  are  finite  in  that  case.  Taking  expectations  in 
Eq.  (40)  with  respect  to  this  distribution  we  obtain  the  following: 


£[*?'(» + 1)]  =  spr/(»)l + ElM  (" + 1)]  -  f'KW  >  o]  = 


=  E[Xl  (n)]  4-  AiE[  length  of  (n  +  l)tk  service  interval  for  queue  i] 


-Pr[X/(n)  >  0] 


(41) 


Now,  in  steady  state,  E\Xi  (n  4-  l)j  =  E[Xj  (n)];  hence 


Pr[X{  >  0]  =  A<E[length  of  service  interval  of  queue  i  in  5,<] 


Since  for  A »  <  Ji\  the  ith  queue  is  stable  in  S3 }  we  have 


1  >  Pr\Xi  >  0]  =  XiPiXi) 


(42) 
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where  /*(A0  denotes  the  expected  length  of  the  service  interval  of  queue  i  in  S3.  Obviously 
this  service  time  is  also  a  function  of  the  other  arrival  rates  A >  /,  but  not  of  the  ones  for 
k  <  y,  since  arrivals  to  the  first  j  queues  have  no  bearing  on  how  often  there  are  transmission 
attempts  from  the  group  of  the  first  j  queues  in  S3  (they  occur  with  probability  1 — ITi=i  (A — P*) 
in  every  slot).  In  any  event,  the  function  f3  is  non-decreasing  in  A,-;  this  is  quite  obvious  since 
more  frequent  arrivals  will  create  more  frequent  collisions  and  less  frequent  successes  by  all 
queues,  including  the  itk  one.  Thus,  for  A,-  < 

/'(AO  <  /'(£')  (43) 

Using  inequality  (43)  in  Eq.  (42)  we  obtain 

1  >  #>  (#) 
or 

/'(«')  <  ^ 

Hence 


Pr[X’  >  0]  =  Ai/'M  <  XifUM)  <  b  (44) 

Now  we  argue  that 

Pr[X{  >  0]  =  Pr[Q{  >  0]  (45) 

Recall  that  Q\  is  the  steady  state  size  of  queue  i  in  S3  at  any  arbitrary  instant  in 
time,  whereas  X(  is  the  steady  state  queue  rise  at  any  arbitrary  departure  point  in  time. 
They  could  conceivably  be  different  random  variables.  However  since  the  joint  queue  size 

[Qi  (n) - Qn(n)]  i®  an  irreducible  aperiodic  Markov  Chain,  ergodicity  implies  that  eq.  (45) 

is  valid. 
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Finally,  by  virtue  of  the  dominance  of  5J  by  £,+1  we  conclude  that 

M  1  <  y  <  M  -  1 

Consequently 

P'lQt  > 0]  <  ^  3+1  <i<M 

This  completes  the  proof. 

Now  we  state  the  main  theorem  of  this  section. 

Theorem  4.1 

The  system  S  is  stable  if 


Xj  <  1  <  3  <  M 


where 


=p;II(1“*)+  pA1-ir)pi  nu-p*) 

i=3+l  ” 


(46) 


<=1 
.  i*3 


k-l 

k*i,i 


and 


bM  =  Mm  =  Pm  II  C1  “  *) 

»=i 

Proof 

By  lemma  4.2  we  know  that  if  A»  <  we  have 


PrlQj  =  0]>l 


i  >y 


(47) 


Substituting  this  in  Eq  (39)  we  obtain 
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Note  that  the  fcj’s,  as  defined  in  Eqs  (46),  (47),  are  less  than  or  equal  to  the  Therefore 
substituting  the  6y’s  in  the  RHS  of  inequality  (48)  we  strengthen  the  bound.  Thus  if  A j  <  bj, 

it  is  also  true  that  A  j  <  /ij  and,  hence,  the  queue  is  stable  in  S'  and,  by  the  dominance 
property,  in  S  as  well.  This  completes  the  proof. 

Notes: 

1.  This  is  an  improved  bound  as  easily  seen  from  the  comparison  of  b{  to  the  value  of 
the  RHS  of  Eq.  (24).  In  fact,  it  is  easily  verified  that  >  p*4. 

2.  The  ordering  of  the  terminals  has  been  arbitrary,  but  the  bound  depends  on  that 
ordering.  The  union  of  the  bounds  of  the  stability  regions  over  all  possible  orderings 
yields  an  improved  bound. 

5.  The  Connection  to  Information  Theory 
The  results  obtained  so  far  for  the  stability  of  the  slotted- ALOHA  system  present  a  striking 
similarity  to  those  concerning  the  achievalbe  region  of  the  no-feedback  collision  channel  [13,14]. 
First  of  all  let  us  refer  to  the  problem  we  studied  here  as  the  ‘‘stability”  problem  and  to  that 
considered  in  [13,14]  as  the  “capacity”  problem.  Let  us  discuss  their  similarities  and  differences: 

“Stability”  problem  -  Determine  the  regions  of  values  of  the  arrival  rates  A*  of  the  Bernoulli 
arrival  processes  for  which  the  discrete-time,  slotted-ALOHA  system  of  M  buffered  users  is 
stable  in  the  sense  of  ergodicity  of  the  underlying  chain,  or  finiteness  of  the  average  queue  size 
(the  two  are  equivalent  here);  that  is,  find  the  values  of  A,-  for  which  3  values  of  p,-,  such  that 
the  system  is  stable.  There  is  no  concern  about  Shannon  capacities,  and  feedback  is  assumed 
to  be  provided  to  the  terminals.  " 

“Capacity”  problem  -  The  object  of  study  is  the  same;  however,  a  different  question  is  asked 
about  it,  namely  at  what  rates  Ri  can  the  terminals  simultaneously  transmit  to  the  common 
receiver  reliably  (in  the  Shannon  sense).  There  is  no  concern  about  stability  of  queues;  in 
fact  the  queues  may  very  well  be  assumed  to  be  infinite  in  order  to  provide  an  inexhaustible 
supply  of  packets,  when  needed,  for  transmission.  There  is  no  feedback,  and  the  terminals  are 
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cot  constrained  to  use  the  ALOHA  protocol  for  transmission,  but  only  some  (any)  protocol 
sequence  of  packet  transmission  and  retransmission. 

To  summarize: 

Both  problems  consider  the  same  channel  (the  collision  channel),  that  requires  a  single 
packet  transmission  to  ensure  successful  reception.  However  they  differ  in  terms  of  the  question 
of  interest  (queue  stability  vs.  reliable  transmission  rate),  and  in  terms  of  the  assumptions  on 
feedback,  arrivals,  and  transmission  protocol. 

Interestingly  the  results  are  remarkably  similar  for  the  two  problems.  In  fact, 

1)  for  M  =  oo  (and  =  Aj  =  ...)  the  “stability"  problem  requires  that  lim  u~m  MXi  < 

a4-.o 

t~l  and  the  “capacity"  problem  states  that  lim*/-..  MR+  <  t~l  is  achievable. 

irf-0 

2)  for  M  =  2  the  “stability”  problem,  as  shown  here  and  in  [8],  requires  that 

\[\x  +  —  I 

and  the  “capacity”  problem  as  shown  in  [13,14]  states  that 

\/Ri  "h  sfR*  £  1 

defines  an  achievable  region. 

3)  for  2  <  M  <  oo  the  “stability”  problem  leads  to  the  bounds  derived  here  and 
elsewhere  [6,7].  These  bounds  can  be  shown  to  be  consistent  with  the  results  of  the 
“capacity”  problem,  which  state  that 


Ri  =  *II(1-*v) 

jV* 

# 

* 

M 

I >•  = 1 

«=1 

defines  an  achievable  region  (capacity  region,  actually,  under  the  constraint  of  using 
a  transmission  protocol  sequence);  namely  that  region  contains  the  region  of  inner 
bounds  obtained  for  the  stability  problem  (as  can  be  shown  quite  easily). 
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These  observations  lead  to  the  following  two  questions. 

1)  Is  there  a  non-superficial  connection  between  the  informationtheoretic  and  the  queue¬ 
ing  theoretic  versions  of  the  problem,  or  is  the  identity  of  the  results  a  mere  coinci¬ 
dence? 

2)  In  either  case,  is  the  conjecture  that  for  2  <  M  <  oo  the  stability  region  is  given  by 

Ai  =  Pi  II*1  “  Pi) 

J*i 

M 

52pi=i 

«=i 


true? 

Whatever  the  answers  to  these  questions  we  may  close  with  two,  hopefully,  useful  remarks. 

i)  Feedback  is  known  to,  in  general,  enlarge  the  capacity  region  of  a  multi-user  channel; 
the  results  seem  to  imply  that,  at  least  for  M=2,  the  ALOHA  transmission  protocol 
makes  poor  use  of  this  feedback  in  that  the  achievable  region  is  not  enlarged,  but, 
rather,  stays  puzzlingly  the  same. 

ii)  The  meaning  of  the  requirement  that  Pi  =  1  in  the  parametric  description  of 
the  capacity  region  remains  unclear,  in  connection  with  the  role  of  the  pi s  in  the 
transmission  protocol  of  the  “stability”  problem.  The  intriguing  nature  of  the  sum 
of  the  values  of  the  p,’s  was  first  noticed  by  Szpankowski  in  [6,7].  Here  we  simply 
note  that  for  the  “stability”  problem  one  can  easily  show  that  if 

M 

^2  Pi  <  1 

*=1 

all  users  can  improve  their  delay  (or  queue-size)  performance  if  they  uniformly  in¬ 
crease  their  respective  p.’s  by  an  amount  given  by 


Ap;  =  P»(l  ~  Pi)  *  At 


(49) 
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where  At  is  a  common  auxiliary  variable,  until  the  sum  hits  the  value  1.  Specifically 
one  can  show  that  any  increment  in  the  values  of  the  p»’ s  according  to  Eq.  (49) 
produces  a  system  dominated  by  the  previous  one  in  the  dominance  sense  that  was 
used  throughout  this  paper.  We  do  not  know  how  to  show  the  improvement  by 
decreasing  the  p,’s  when  Y^iLi  P »  >  *• 

6.  Conclusions 

In  this  paper  we  have  used  a  new  technique  involving  stochastic  dominance  to  study  the 
ergodicity  region  of  the  discrete-time  slotted- ALOHA  system  with  a  finite  number  of  buffered 
users.  We  have  obtained  the  region  exactly  for  M=2  and  derived  inner  bounds  to  it  for 
M  >  2,  that  improve  upon  earlier  ones.  The  results  look  remarkably  similar  to  those  of 
the  no-feedback  collision  channel  and  thus  lead  to  speculation  about  a  deeper  relationship 
between  the  queueing-theoretic  and  information-theoretic  approaches  to  the  problem.  They 
also  suggest  a  conjecture  about  the  precise  boundaries  of  the  stability  region  for  M  >  2.  The 
dominance  approach  introduced  in  this  paper  represents  a  natural  and  simple  concept  that 
may  be  of  further  value  to  the  study  of  other  related  problems. 
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Appendix  1 

In  this  appendix  we  establish  the  dynamics  of  the  ith  queue  in  system  S’,  as  defined  in 
the  general  formulation  of  section  4,  and  obtain  the  stability  criterion  that  is  used  repeatedly 
in  all  sections  of  the  paper.  Specifically  the  objective  is  to  show  that  in  S’,  if  A»  <  /*?’,  the  ith 
queue  is  stable  in  the  sense  that 


P’,o  =  P'lQi  =  0]  >  o. 

Consider  the  ith  queue  of  S’  at  time  n  and  the  system  equations  for  the  queue-size 
probabilities.  These  are 

P/,o(n)  =  ~  1)  +  A*M-tl(n  -  l)P/’x(n  -  1)  (Al.l) 

P/.i(n)  =  (n  -  !)  +  ~  l)]p/’i(*»  -  1) 

+A<P/;0(»  -  1) 

+A»<1 (»  -  l)p/’2(»  “  0  (A1.2) 

Pi,kin)  =  ~  ^)Pi,k-li.n  -  1) 

+*»•£<,* (n  -  J)  +  A‘<t(n  -  1  )Plk{n  -  1) 

+*»<*+ 1(»  “  l)pi,k+ 1(»  -  1)  (A1.3) 

Obviously  the  “death-rates”  p£fc(»)  are  functions  of  the  status  of  the  other  queues,  which  is 
why  the  queues  are  coupled  and  difficult  to  anlayze.  However,  the  above  equations  are  valid 
conditioned  on  the  status  of  the  other  queues.  There  is,  in  our  notation,  an  explicit  dependence 
shown  on  the  size  of  the  ith  queue.  Thus 


4fc(n)  =  Pr[P?(n)  =  l|g?(n)  =  fc] 


Suppose  now  that  P/>Jk(n)  and  Mi,k(n)  posseaa  limits  as  »  -♦  oo;  that  they  indeed  do  will 
be  established  a  little  later.  By  considering  the  limits  of  the  system  equations  we  obtain 
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(A1.4) 


PL  = 

+*&7i  +  *it*i,iph 

+^«Mf,a^a 


(A1.5) 


PU  =  XKk-iPlk-i 
+X*fii,k+lPitk+l 


(A1.6) 


These  equations  can  be  solved  to  yield 


pi  —  pi  X* 


(j41.7) 


•  A*  1  TT  t 

P^  —  pi  ___  TT  iV 


(A1.8) 


and  hence 


i  +  tV+  £  A-  n‘(¥r) 

*<**•,1  k= 2  *=l  **<.**» 


(Al.9) 


Based  on  Eq.  (A1.9)  we  conclude  that: 
1.  Queue  i  is  stable  in  the  sense  that 


Pr[Q?'  =  0]  =  /V.o>° 
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(>11.10) 


A* <  j£™X* 

*■^00 

2.  Queue  i  b  unstable  in  the  sense  that 


if 


Pr[Q\  =  oo]  =  =  1 

Xi  >  lim  ftk 
*-♦00  * 


(Al.ll) 


Note  also  that  Eqs.  (A1.8)  and  (Al.9)  imply  that  under  the  stability  condition  the  moments 
of  the  queue  sues  exist  as  well. 

We  now  wbh  to  show  that  Xi  <  ft  implies  stability.  If  <  ft,  one  of  the  following  b 

true: 

Case  1  -  A*  <  £ {  <  lim  ft  k 
Case  2  -  A<  <  lim  ft  .  <  ft 
Case  3  -  Ai  =  lim  ft  .  <  ft 
Case  4  -  JUm^^  <  A<  <  ft 

Cases  1  and  2  satbfy  Eq.  (10)  and  hence  the  queue  b  stable. 

Case  4  cannot  occur,  because  if  it  does  then  by  Eq.  (Al.ll)  lim  P?k  =  1,  which  imples 

k— *co  *’* 

that 


ft  =  ^  ft,k^i,k  =  ft,* 


k= 0 


which  contradicts  the  inequality 


lim  ft  k  <  ft. 

k—co  *>*  ’ 
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